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A dissipation rate, which grows faster than any power of the wave number in Fourier space, may be scaled 
to lead a hydrodynamic system actually or potentially converge to its Galerkin truncation. Actual convergence 
we name for the asymptotic truncation at a finite wavenumber kc above which modes have no dynamics; and, 
we define potential convergence for the truncation at kc which, however, grows without bound. Both types 
of convergence can be obtained with the dissipation rate ii[cosh(k / kc) — 1] who behaves as k^ (newtonian) 
and expl/c/Zcc} for small and large k/kc respectively. Competition physics of cascade, thermalization and 
dissipation are discussed with numerical Navier-Stokes turbulence, emphasizing on the intermittency growth. 
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Energy cascade in three-dimensional (3D) turbulence is ac- 
complished by hierarchically generating smaller and weaker 
eddies (which are eventually eliminated by dissipation,) while 
thermalization in a conserved system is by sharing (general- 
ized) energies among all excitations equally. Fermi et al. [IJ 
found some nonlinear systems did not simply thermalize as 
intuitively expected. However, it has been recently shown 
[2 1 that even some dissipative systems do partially thermalize, 
somewhat counterintuitively again. The idea is simply that the 
normal dissipation is just an incomplete truncation. The trun- 
cated Euler equations are a Liouville system and may thermal- 
ize to an equipartitioned k'^ spectrum in the three dimensional 
nonhelical case El 01. It is then proposed to explain the flatter 
spectrum in between the inertial and dissipation ranges in fluid 
turbulence, usually called a bottleneck (see L2J and references 
therein,) as partial thermalization. 

The notion of partial thermalization was proposed with the 
use of a high power a of the Laplacian (hyperviscosity) in 
the dissipative term of hydrodynamical equations. More pre- 
cisely, consider the hydrodynamic + /i(^)^"]u(k) = 

B(u, u) for the velocity field u in a cyclic box represented in 
Fourier (wave number k) space, with kc being off-lattice, and 
B(-, •) is the nonlinear term. With a ^ oo while the positive 
but finite ji and kc being held fixed, the dynamics corresponds 
to the Galerkin-truncated-at-Zcc equation dtv = n/e^5(v, v)\ 
The truncation projector like is defined as a low-pass filtering 
operator which keeps the harmonics with wavenumber less 
than kc and sets the other ones to zero; v — Vl]^^u. We thus 
see that equilibrium statistical physics plays a role in compe- 
tition with the unique dynamics controlled by the structure of 
dynamical euqations. On the other hand, inertial range of tur- 
bulence, characterized by a constant energy flux, is nonequi- 
librium, and it has been observed that some aspects of the 
intermittent cascade physics of turbulence also work beyond 
the inertial state (see, e.g., O.) The interesting lesson we thus 
have learned is that turbulence is so generous as to allow the 
persistence of both these two extreme cases of statistical me- 
chanics ! 

As, from the above, dissipation is responsible for both 
cascade and thermalization, what then is the essence of the 



dissipation mechanism needed to obtain convergence of the 
dynamics to a Galerkin truncation, and that thermalization? 
What if the asymptotic truncation wavenumber grows without 
limit, in some particular way? And, how does thermalization 
affect the dyanamics of the flow structures which are gener- 
ally counted for turbulence intermittency? Answers to them 
are essential to the general understanding of nonlinear dissi- 
pative systems and need relevant details of explicit examples. 
We start with a brief review discussion with our results: 

a) We notice from Frisch et al. Q that the dissipation 
rate needs to grow faster than any power of k to converge 
to a Galerkin-truncation operator of the system. Exponential 
growth, among others, of the dissipation rate will be shown to 
be also able to lead to such Galerkin truncation of a hydrody- 
namic system. Dissipation rate 

V(k) = /j.[cosh{k/kc) - 1], (1) 

now called coshcosity, which grows exponentially for large 
k is expected | lOJ to be enough to tame the solutions to be 
not only analytic but also entire, is a perfect model for our 
purposes here. 

b) It was argued that the Galerkin-truncation for a hydro- 
dynamic system has a mathematical consistency problem in 
the sense that the symplectic structure and Casimirs are bro- 
ken, and a self-consistent sine truncation (for two-dimension 
hydrodynamics) is then proposed |8|. One should notice that 
with the presence of growing dissipation rate at high modes 
the Galerkin truncation is actually physically relevant O. 
And, we do not even know what the Euler equations are ||4j|9l 
so that discussing convergence to 3D Euler equations seems 
pre-mature. However, we may numerically obtain some indi- 
cation of what would happen when kc goes to infinity in some 
specific way. We denote one type of asymptotic behavior as 
"potential convergence" to Galerkin truncation. 

c) Boyd | 6 | studied various dissipation models to observe 
their effects on a shock structure. The statistical consequences 
on turbulence should then also be investigated in detail. In tur- 
bulent flows, coherent flow structures 1 7 1 appear to be "spotty" 
as the insight of intermittency given by Onsager (H. We will 
present numerical Navier-Stokes results showing that small- 
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scale coherence is destructed to be mistily and uniformly dis- 
tributed with a reduction of intermittency. 

The local dissipativity r{k) = in Eq. |lj) cor- 

responds to a in hyperviscosity V{k) = When 
a is increased without limit hyperviscous V{k) grows faster 
and faster, approaching zero below kc and infinitely large 
above kc (so that energy can not "tunnel" through any more.) 
Now r{k) from the coshcosity model grows with k with- 
out limit, implying the potential of leading to the Galerkin 
truncation in a unique way. As the dissipation time scale is 
[cosh(/c/ kc) — with the eddy turn-over time fixed or 
varying slowly we can increase the dissipation wave number 
kd by decreasing with fixed kc. This way, kd essentially 
becomes a truncation wavenumber kc- Since kc is unlimited 
with /i ^ 0, we say potential convergence to distinguish from 
the actual convergence with finite kc as in 1 2 1 and in the other 
parameterization of coshcosity immediately below. Similar 
results may be obtained by letting a be some increasing func- 
tion of kc in the hyperviscosity model. Actually results below 
seem to indicate that, if r(/c) grows with k fast enough, a lim- 
iting system with conservative dynamics may be obtained. 

We now study how to obtain actual convergence to the 
Galerkin truncation at some finite wave number kc with the 
coshcosity model. Following the idea of the estimation of the 
dissipation wave number kd in the above, we should have the 
cross wave number kc decrease with /i to prevent kd from 
growing but to make it approach some finite value. This may 
be accomplished with the following observation: The right 
hand side of Eq. ^ is essentially for large k, so by 

letting kc — — /cg/ In /i, as /i ^ 0+, 

,l-fe//CG 



V(k) - /i' 



(2) 



V(ki) then goes to zero for k < kc and to infinity for k > kc, 
which may lead to actual convergence to the Galerkin trunca- 
tion. Noticing that when k is small, the coshcosity reduces to 
normal viscosity by taking the leading order term in the Tay- 
lor expansion, for /i not small enough we need to make some 
subtle adjustments, e.g., to take kc = kc/^i^if^ + — 1) 
with K, = l//i + las used below, to make the peak of the 
energy spectrum approach kc in a consistent way. 

To further illustrate the ideas in the above we first per- 
form the integrations of eddy-damped quasi-normal Marko- 
vian (EDQNM) equation for turbulence energy spectra as in 
121 . The equations read |11 1 



[dt^2V{k)]E{k,t) 



Afc 



dpdq Okpq X 



b{k,p,q)^E{q,t) [k^E{p,t)-p'E{k,t)] ^ F{k) 
pq 
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b{k,p, q) = yixy^ z^) , Okpq ^ ^ , 



fik = V{k) + A 



Jo 



p^E{p,t)dp 



Here, E{k^ t) is the energy spectrum at time t, F{k) is the 
(spatial) forcing spectrum, x, y, z are the cosines of the angles 



of the triangle with sides k, p, q and A is a constant related 
to the Kolmogorov constant. We use the same forcing, dis- 
cretization, and, mixed time-marching and iteration schemes, 
except that the hyperviscosity is now replaced by coshcosity. 
Two types of parameterizations introduced above are applied 
resulting in Fig. [T] where the stationary spectra show respec- 
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FIG. 1: EDQNM spectra with coshcosity ([T}: kc = kG/\n{i<i + 
— 1) and K = + 1 (lower blues lines), and kc = constant 
(upper red lines.) Fixed values of kc and kc are 10^ and 10^ re- 
spectively. The line named viscous is the same spectrum in | 2 1 with 
hyperviscosity V(k) = for a = 1 case (actually normal 

viscosity now) with fi — 1 and kc — 10^. 



tively, as /i ^ 0, actual convergence to Galerkin-truncation at 
a finite wave number kc (upper red lines) and potential con- 
vergence to truncation at — In /i (lower blue lines.) We see 
clearly a inertial range followed by a little bit of sec- 

ondary bottleneck (a largest overshoot of 3% for /i = 10~^^ 
also can be seen at around k = 2000 if zoomed in,) a pseudo- 
dissipation range, a thermalization range and finally a dissi- 
pation range. The traditional assumption of inertial scaling 
going straight down to dissipation scale for estimating dissi- 
pation scale is vitiated by partial thermalization. When the 
thermalization is strong, the dissipation rate changes drasti- 
cally around kd which can be conveniently estimated numeri- 
cally by taking V{kd) ^ 1, as long as the eddy turnover time 
is not extremely far from order one (which is generally the 
case in numerical simulations) and then can be captured by 
V{kd + 5k) with \5k\ being relatively small. If kc is kept 
constant and only ji itself varies, following the above phe- 
nomenology, we have kd ~ — /ccln/i as designated by the 
vertical dotted line (?^ 92103) for /i = 10~^^ case. The sys- 
tem has a dissipation wavenumber kd which goes to infinity as 
fi and become essentially kc- One might conclude from 
the upper lines for potential convergence that a ("directional") 
limit of conservative spectral dynamics may be obtained in 
this way 1 19 |. 

Knowledge of how the thermalization and dissipation de- 
pend on the parameters may be useful for turbulence simula- 
tions and modeling. With fi large, the majority of the spectrum 
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falls in the regime where coshcosity reduces to normal viscos- 
ity. Here we compare the normal viscosity case, exactly the 
same as the a = 1 case in |2|, and a coshcosity case with 
/i = 10^. The latter may be of practical interest: The begin- 
ning of the dissipation range and below is mostly as with a 
normal viscosity, to which the coshcosity is reduced for small 
k. Only at higher wavenumbers the faster growth of the dissi- 
pation rate is felt, with the fluctuations there being low and of 
little practical significance. Such case may be used to avoid 
wasting resolution without developing a serious bottleneck as 
also proposed in 1 10 |. 

Having the concepts of actual and potential convergence 
to Galerkin truncation in mind, together with the knowl- 
edge of EDQNM spectral dynamics, we then proceed to in- 
tegrate the 3D Navier-Stokes dynamics (with the nonlinear 
term B{-, •) = -ikmPij Ep+q=k ^i(p)^rn(q), = -1 and 
the transverse projection operator Pij = {Sijk'^ — kikj)/k'^.) 
We explore the physical contents of cascade, thermalization 
and dissipation concerning the statistics of velocity increment 
ArU = u{x + r) — u{x) which is a deep quantity |4, 9|. It 
suffices for our main purposes here to report the results of the 
smallest /i (= 10~^^) case from a 512^ simulation in a cyclic 
box of period 27r with kc = 120.67. 

Fig. [2] shows the structure functions Sp{r) =< \Aru\P > 




FIG. 2: Structure functions: Bottleneck scales are roughly desig- 
nated by two vertical dashed lines. The analytical scalings for p = 2 
and 14 and the scale-independent scaling (oc r°) are also plotted for 
reference. 5* is a number used to shift the lines properly for the mir- 
rored log-log plot against S3 ■ 



(with < • > meaning statistical average) against scales (left) 
and against the third order structure function (right). We see 
that bottleneck scales in physical space now also clearly show 
up as designated by the corresponding pairs of vertical dashed 
lines. The fast damping beyond a scale r] in physical space 
can be related to the spacial discretization with grid width r]. 
Identical independent Gaussian distribution, as of the equilib- 
rium of Galerkin-truncated system in Fourier space, leads to 
also Gaussian distribution with negligible (decreasing with to- 
tal number of grid points as can be shown with simple algebra) 
non-diagonal elements in the covariance matrix of the spa- 



cially discretized velocity (c.f. the reference scaling on the 
left plot.) Clearly from the right part of the plot, Sp against Ss, 
for the extended self-similarity (ESS - claiming power-law re- 
lations between structure functions Sp{r) ~ [Sq{r)]^P''i il3J ,) 
the bottleneck regime is a transitional range from analytical 
subrange scaling (ASRS, with the data highlighted with blue 
circles,) = p to the inertial subrange scaling (ISRS) (p, 
with the data highlighted with red triangles: Although, from 
the left part of the figure we see that the pure analytical range 
is barely resolved (cf. the analytical scalings shown in the fig- 
ure for p = 2 and 14 for comparison [20 1) and a pure inertial 
range scarcely emerges, the right part of the figure shows that 
ESS works well |21 1 both for analytical scaling and inertial 
scaling but fails in the bottleneck scales - bulges can be found 
here for various p with careful observation. Using ESS, the 

analytical scaling is measured, through Sp = ApS^^ , to be 
= 1.0076p. The inertial scaling exponents (p show no ev- 
idence of deviation from those measured in normal fluids. So, 
in our simulation, some asymptotic statistics in both inertial 
and analytical ranges are already available. 

Incompatibility of ESS in the bottleneck region has actu- 
ally already been found in Ref. tT4\ with the careful obser- 
vation of the intermittency growth. As usual we thus measure 
the intermittency by the flatness factor of velocity increments 
F4(r) = 6'4(r)/[5'2(r)]2. As shown in Fig. ^ an obvious re- 
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FIG. 3: Flatness factor for the velocity increments: The same bot- 
tleneck regime as in Fig. |2]is designated by a pair of dashed lines. 
A surrogate = 3.3871, by ESS, is measured for the asymp- 

totic flatness factor in the analytical range S^j S\ — a^ja^. The 
flatness inertial-range scaling (measured by ESS), with exponent 
C4 — 2^2 = —0.1129, is also shown to fit well to the data of Navier- 
Stokes with normal viscosity ("viscous" - reproduced from 1 14 1). 

duction of intermittency does present in the bottleneck regime, 
in comparison with the tiny "lull" with normal viscosity {cf. 
the "viscous" data reproduced from |[T4ll .) Coherent vortex 
motions as a signification of the self-organization of the flow 
are randomized at such scales. Once r goes into the analytical 
range with Spiv) = apV^, the flatness factor should be a con- 
stant depending on the settings of the flow, say, the Reynolds 
number (see, e.g., L16J ). The much smaller flatness factors 
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in the far dissipation range for our coshcosity data than those 
with normal dissipation as shown in Fig. |3]may be further 
explained as follows. Strong damping of high modes push 
the complex singularities (if any) further away from the real 
axis(es) (even to infinity as proposed in 1 10|) and/or weaken- 
ing them. This then reduces the intermittency by the Frisch- 
Morf argument for the origin of intermittency 1 15 |. 

It is then intriguing to see the effects on the flow config- 
urations of the competition between thermalization and self- 
organization. The top-left panel of Fig. |4] shows the larger- 
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FIG. 4: Direct volume rendering ifTSl of du/dx\ The top-left panel 
is the filtered field keeping only the modes under the thermalization 
range; the top-right panel is the total field; the lower panels present 
the linear-log plots of the histograms (of the corresponding fields 
of the upper panels) and the transfer functions. The du/dx val- 
ues where the transfer functions (graphed by the lines connecting 
the opacity-control points) for opacity beginning to be nonzero (for 
visibility) is approximately the values where contributions to flatness 
factors peak. Changing the transfer function will make the picture 
look more clear or cloudy instead of the cleanly spotty versus mistily 
uniform properties in the two renderings. 



scale structures are "'spotty' distribution of regions in which 
the velocity varies rapidly between neighboring points" as un- 
derstood by Onsager IH. Because, due to the eddy viscosity 
caused by strong thermalization of high modes which are fil- 
tered out now, lower modes below the thermalized range work 
like a normal fluid 1 17 |. These structures are embedded in the 
very-small-scale dissipative structures which are almost mist- 
ily uniform as shown in the top-right panel for the total field 
(including those thermalized modes.) 

We thus have shown by applications of the coshcosity 
model that a dissipation rate growing faster than any power 
may be parameterized to be of either potential or actual 
Galerkin truncation convergence capabilities. The dissipation 
wave number kd in either case approaches to the Galerkin- 
truncation wave number kc which grows without bound and 
is a finite value respectively. The EDQNM calculations for 
the latter case may indicate a limit of conservative spectral dy- 



namics. Measurements and analyses of intermittency growth, 
ESS properties and careful visualization of the flow fields 
from a direct numerical simulation of Navier-Stokes equa- 
tions detail the following understanding: Due to dissipation at 
small scales, turbulent flows are composed of competitive and 
complementary cascade and thermalization, which is pictured 
as self-organization and destruction of coherent structures at 
successive scales. 
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